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Abstract: We compute the short-circuit diffusion current of excitons in 
an organic solar cell, with special emphasis on fluorescence losses. The 
exciton diffusion length is not uniform but varies with its position within 
the device, even with moderate fluorescence quantum efficiency. With large 
quantum efficiencies, the rate of fluorescence can be strongly reduced with 
proper choices of the geometrical and dielectric parameters. In this way, the 
diffusion length can be increased and the device performance significantly 
improved. 
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1. Introduction 

A fundamental obstacle to the efficient conversion of solar radiation into electric current by 
organic material is the short diffusion length of photo-generated excitons. These need to survive 
long enough to reach a dissociation site, where they can separate into holes and electrons. 
Hence the thickness of the donor and acceptor layers in bilayer hetero-junctions is severely 
limited. That problem is overcome in bulk hetero-junctions, but it is then replaced by the fact 
that electrons and holes need to percolate through an intricate structure. They are thus less 
mobile and can be trapped in dead ends (T). The aim of this paper is to show that the exciton 
diffusion length can significantly be modified by the micro-cavity formed by the various layers 
in an organic solar cell, and that a proper understanding of this effect could lead to a better 
performance of bilayer hetero-junctions. 

Most exciton transport models boil down to the diffusion equation 

= L 2 ^--p+g(z), L 2 =Dt, (1) 

in the photoactive material, where p is the exciton concentration, D is the diffusion constant, T 
is the lifetime, and g(z) is the source term. Different modeling considerations lead to different 
forms of g(z) and different boundary conditions ||2j-j4] . A comparison between the main models 
concluded that they all qualitatively reproduce the photovoltaic response of a given organic 
solar cell but that none was fully accurate Q. More recently, (Q~|) was used to experimentally 
determine the diffusion length L from the optical response of solar cells J^[8]. 
However, if fluorescence is taken into account, ([TJ should be modified as 

= L 2 ^--b(z)p+g(z), L 2 =Dt, (2) 

where the factor b(z) is the decay rate normalized to the bulk value and is now space-dependent. 
Thus, the diffusion length is locally corrected by the factor 1 / y/b{z). The function b{z) is de- 
termined by the optical thicknesses of the layers that make up the solar cell and the wavelength 
of the fluorescent light (9j. It can differ significantly from unity. The fraction of the exciton 
decay rate that is due to fluorescence, or 'fluorescence quantum yield' is usually assumed small 
in organic solar cells. In this paper, we specifically address the case of a large quantum yield 
and explore the possibilities afforded by this process in improving device performance. 

That fluorescence is affected by the device geometry is already well accounted for in the 
design of organic light-emitting diodes (OLEDs), see for instance |[T0UT8ll . For OLEDs, the 
diffusive transport of excitons is less of an issue than for organic photovoltaics, although this 
aspect has been considered in H 1 911201 . Indeed, it is obviously desired that OLED excitons fluo- 
resce and they are expected to do so near the interface between donor and acceptor molecules. 
For photovoltaic devices, we have recently drawn attention to the usefulness of controlling flu- 
orescence in Schottky solar cells lETI . To this end, we had computed fluorescence rates using 
Fermi's golden rule, and had therefore neglected dissipative processes, such as non radiative en- 
ergy transfer to metal electrodes. In this paper, we model excitons as classical dipole emitters, 
allowing us to correctly treat dissipation in the electrodes. The radiation of an electromagnetic 
dipole in a stratified medium has been the subject of many theoretical works ||9]|22 ti27l . all 
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Fig. 1. Normalized fluorescent decay rate as a function of distance from a thick Ag electrode 
in a uniform medium of refractive index n = 1.3. Thin lines: ideal, lossless electrode, was = 
6.37/. Thick line: real electrode, «Ae = 0.04 + 6.37/ 1341 . Vaccuum emission wavelength: 
900 nm. bu : parallel excitons. perpendicular excitons; 

building from Sommerfeld's analysis of a radio-wave antenna above the Earth [28]. Finally, let 
us note that the idea of fluorescence management in photovoltaic cells was already put forward 
in ||29l , but it was discussed only at a qualitative level. 

One may approach fluorescence management from two limit situations: (i) an emitter inside 
a perfect, sub-wavelength, cavity and (ii) an emitter in front of a single mirror. From the former, 
we know l30l that complete suppression of fluorescence is possible. However, no exterior radi- 
ation is admitted inside such a cavity. The latter configuration, being open to the environment, 
appears to be more relevant to the present discussion. We shall call it the 'Drexhage config- 
uration', in reference to Drexhage's famous study of fluorescence near a mirror H31JI321 . The 
function b(z) can be obtained from the power radiated by an electromagnetic dipole in this con- 
figuration [28 33 1 . To this end, one should distinguish excitons that are oriented parallel to the 
electrode (b\\) from those which are perpendicular to it (b±), see Fig.Q] 

Let us assume that the photoactive material has a refractive index n = 1 .3 and that the exciton 
radiates at a vacuum wavelength of 900 nm. At that wavelength, the refraction index of silver 
is nAg = 0.04 + 6.37/ 134| . If silver were lossless, i.e., if «a s were purely imaginary, then it is 
evident from Fig. Q]that interesting prospects of device improvements would be afforded by 
the low values of bn(z) near the electrode, z = 0, and of b±(z) near z = 250 nm. Note that, 
since the functions b\\ (z) and b±(z) vary in opposite ways, it appears that taking advantage of 
space-dependent fluorescence requires one to control exciton orientation. Let us next look at the 
decay rates in the vicinity of real silver. Due to dissipation, the decay rates diverge as the exciton 
comes close to the electrode. On the one hand, b\\ (z) still has a well-pronounced minimum a 
few tens of nm away from the electrode. The lower the value of n, the lower the minimum, 
and with n = 1 .3, the minimum of b\\ (z) ~ 0.25, representing a doubling of the diffusion length 
in (f2|i. On the other hand, Fig. Q] shows that perpendicular excitons are affected by dissipation 
over a significantly longer distance than parallel ones. In this particular configuration, b±(z) is 
nowhere brought close to zero; micro-cavity effects from more complicated layered geometries 
are necessary to achieve this. 

Returning to parallel excitons in Fig.Q] we note that the minimum of (z) is very close to the 
silver electrode. On the other hand, g(z), being proportional to sunlight intensity, is also close to 
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Fig. 2. Solar cell geometry and coordinate system. HBL: hole-blocking layer. A: Acceptor. 
D: Donor; EBL: electron-blocking layer. 

zero near the silver electrode due to interference. The advantage of having a small b(z) is thus 
lost in that situation. Hence, in the search of an optimal design, both functions b(z) and g(z) 
must be carefully monitored. This is in line with previous work demonstrating the sensitivity of 
g(z), and hence device performance, on the geometrical parameters 13511361 . 

At first sight, it may seem that decreasing the spontaneous emission in a solar cell device 
would also affect stimulated absorption. Indeed, Einstein's coefficients of spontaneous emission 
and stimulated absorption, respectively A and B, are related in free space by A = (87r/iv 3 /c 3 )B, 
where h is Planck's constant. In fact, the relation between theses two constants is more generally 
given by A = M(v)hvB, where M(v) is the effective spectral mode density per unit volume. 
This mode density can be considerably less in a confined environment than in open space, 
allowing one to reduce spontaneous emission while leaving absorption undisturbed. In the same 
vein, one should note that, in a solar cell, stimulated absorption is primarily caused by photons 
coming at normal incidence from the sun, while spontaneous emission occurs in all directions. 
In the latter process, the details of the resulting electromagnetic field, and hence, the coupling 
of the emitter to that field, depends on the direction of emission as well as the position within 
the microcavity. This also affects M(v). 

The rest of the paper is organised as follows. In section|2] we discuss the mathematical model 
© in more detail. Section[3]presents our numerical results. In § |3.1| we consider a configuration 
that is favourable to parallel exciton transport, while in § |3.2| we present a configuration that is 
suited to perpendicular exci tons. Finally, in section|4] we present our conclusions. 

2. Model 

2.1. Diffusion 

A typical bilayer organic solar cell architecture is sketched in Fig. [2] The photoactive materials 
occupy the regions z-i < z < zo and zo < Z < Z\- In each of these, the density of excitons, p' 
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where Z), is the diffusion constant, T, is the bulk value of the exciton lifetime, a,- is the sunlight 
absorption coefficient, 0,- the quantum efficiency of exciton generation, N is the incoming pho- 
ton flux, and g(z) is the distribution of sunlight intensity inside the device. g(z) is computed 
using transfer matrices and assuming a normally incident wave with unit photon flux at a spe- 
cific wavelength, A s . Given the complex index of refraction «, . s in region i at X s , the absorption 
coefficient a, is given by 4nlm(n^ s ) / X s . 

In each photoactive region, we normalize p' as 



p' = a^NXiP, 



(4) 



giving 



0: 



f(z), 



Lf=D,T„ 
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where p, and g are dimensionless. We assume that complete exciton dissociation into elec- 
trons and holes occurs at the Donor/Acceptor interface. The hole-blocking layer (HBL) is as- 
sumed to block exciton current but to let electrons flow easily. Similarly, electron-blocking 
layers (EBL) are supposed to let holes through but to block exciton current. Hence, (0 should 
be solved subjected to the boundary conditions 
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(6) 



In the absence of blocking layers, the electromagnetic losses in the electrode make b(z) diverge 
as z Z±i, resulting in a vanishingly small exciton density at these locations, independently of 
the imposed boundary condition. 

Let us note that more detailed transport models than mere exciton diffusion exist, particularly 
drift-diffusion models 1071 which include electrons, holes and static electric field distributions 
within the cell. However, exciton transport asymptotically decouples from the rest of the system 
and governs the device dynamics in the limit of small exciton mobility ll38l . 

If all excitons reaching the interface at zo disintegrate into a pair of electron and hole and 
all of these reach the electrodes, then the short-circuit current is the sum of the two diffusion 
currents at zo- Using |@}, it can be expressed as 



dp 
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A being the area of the interface. We may group the terms in the two diffusion currents above 
as 



I i = ±(AN)a i <S> i L 2 i ^- 
dz 



(8) 



In this expression, AN is the number of photons falling down on the device per unit time, 01,0, 
is the fraction of them per unit length which are converted into excitons, and the diffusion con- 
stants and bulk lifetimes only appear through the diffusion length L;. Let us note that 7 SC / (AN) 
is the external quantum efficiency (EQE) under the short-circuit condition. We thus have 



EQE S( 



i T 2 dp 
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In this work, we only monitor the short-circuit current as a figure of merit for a given configura- 
tion. Indeed, this can be directly obtained from the solution of ©, without additional modeling 
assumptions. 

Let us assume that 0, ■ = 1 . We can obtain a useful estimation of EQE SC if the photoactive 
region is sufficiently thin that we can replace b(z) and g(z) by their average values in (0: 

OttLf^--bp+g, b = -J b(z)<k, g = \( g(z)dz, (10) 



' dz 2 ' d Ja+d ~ ' d Ja+d 

where d is the total thickness of the photoactive region. If the Donor and Acceptor material have 
the same diffusion length L, the maximum short-circuit current is obtained when both materials 
have the same thickness d /2. A straightforward calculation then shows that 

2L / dVb\ 

EQE SC w agd x — ^ tanh I I = r\ A x r\ D d 1) 

where t\a is the fraction of the photon flux that is absorbed by the photoactive region and rjo is 
a diffusion efficiency, in the same spirit as in [TJ . 

2.2. Decay rates 

We now turn to the computation of the normalized decay rate b(z). Its expression was deter- 
mined for a dipole within an arbitrary stack of parallel layers of isotropic materials by Chance, 
Prock, and Silbey [ 9 1 . The case of a dipole embedded in a uniaxial medium sandwiched between 
two other uniaxial media with aligned extraordinary axes was treated in J25). This last case can 
easily be generalized to an arbitrary stacking of aligned uniaxial media, as we do below. Note 
that the Green dyadic was computed for a uniaxial multilayered media in (23), although with no 
particular emphasis on b(z). For the sake of completeness, we provide an alternative derivation 
of b(z) to that given in ||25l and give the relevant components of the Green dyadic in appendix 
to this paper. 

When computing the fluorescence in the photoactive layers, it is necessary to assume that 
the refraction index is purely real in these particular layers. Otherwise, the power required from 
an electric dipole to sustain harmonic oscillations is infinite and the validity of the classical 
formulas given below is uncertain. Fortunately, a frequent feature of organic materials is that 
the exciton energy level is well below the LUMO, with a binding energy ranging between 
0.1 eV and 1 eV. Hence, the exciton radiation is only weakly absorbed by the surrounding 
material and the refraction index at the exciton frequency is mostly real. It appears reasonable, 
then, to neglect the imaginary part of the refractive index in the donor and acceptor material at 
the exciton frequency. 

Let us assume that the layers in the device are uniaxial and that they all have their extraordi- 
nary axis in the z direction. The relative permittivity tensor reads 
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Given the vacuum wavenumber &o, the dispersion relation for extraordinary waves, or /7-waves, 
is 

^ -fr*f' (13) 

between the z-component, k e -j, of the wave vector and its projection hi on the (x,y) plane. On 
the other hand, ordinary waves, or s-waves, satisfy the dispersion relation 



k ^ = Je^kl-k\. (14) 



With these notations in mind, the normalized decay rate for an exciton oriented in the z direc- 
tion, and located between zo and z\, is given by (see l25ll or the Appendix) 




where q is the fluorescence quantum yield, 

R~ p = R p exp[2ik e , zA {z-zo)} , R'l = R';exp[2ik e -A ( Zi -z)} , (16) 

and Rq and R p are the coefficients of reflection from the layers below and above the exciton, 
respectively. 

To compute the reflection coefficient R 1 ^ of a /j-wave incident from a medium with permittiv- 
ity tensors £, onto a multilayered medium with permittivity tensors e, + 1 , £, + 2 ,...,£#, £#+ 1 and 
thicknesses dj+i, c/,+2, . . . ,djv we use the downwards recurrence 



j ~ l 1 + fiJL ! j exp {2ik e ,z,jdj) ' 



where Rf, is the p-wave reflection coefficients between half-spaces filled with media ; and j: 

nP _ ke,z,j£j,x ~ k e ,z,j£i.x 

'j L r ■ L r ' 

For an exciton with dipole moment in the (x,y) plane, and z between zq and z\, the decay rate 
is given by the formula 

oil (z) = 1 + 



4 3£ La - + £ L - 



xRe -/ T 9 V *„L +1 / Wl • (20) 



Here, 



fl = /$exp [2ik , Z;1 (z - zo)] , R\ = R\ exp [2ik <>jZ ,i (zi - z)] , (21) 

where Rq and R\ are computed by a similar recurrence as for /5-waves: 



R) Kj • A-expiI/Ay,, _ kp £l -ko M 
R N -R W , R ^- 1+RUjR s exp( 2 ikozjdj y R iJ~k~ Tk— ■ 



The computation of £>x and ftp in the region z~\ < z < zo follows the same pattern. 
Finally, for randomly oriented excitons, b = + ^b±. 

3. Numerical results 

We solved Eq. © in the two active regions depicted in Fig.|2]for various device configurations. 
The bottom electrode was assumed to be Ag in all cases and the transparent electrode was taken 
to be either ITO or a thin Ag layer. Note that other metals, such as Au or Al, could equally 



be used in the simulations; from an electromagnetic point of view, they influence the device 
performance in the same way as Ag. The refraction indexes of Ag and ITO were taken from 041 
and 11391 . respectively. The rate of production of excitons, g(z), was computed at a wavelength 
A s = 750 nm, where the solar photon flux is highest. On the other hand, fluorescence was 
computed at a wavelength A = 900 nm, corresponding to an exciton binding energy of 0.275 e V. 
For each simulation, the layer thicknesses must be specified as well as their refractive indexes 
at A and A. s . For the exciton- and hole-blocking layers, we assumed identical real refraction 
indexes at A and A. v , ranging between 1.3 and 2.8. In the Donor and Acceptor materials, an 
imaginary part was added to the refraction indexes at A. s in order to account for a prescribed 
sunlight absorption coefficient. In all cases, we assumed an absorption length, a -1 , of 70 nm. 
In addition, the quantum yield and the diffusion length should be prescribed in each photo- 
active layer. Having so many parameters to vary, we imposed e x = E- everywhere despite the 
anisotropy implied by assuming either b = b\\ or b = b±. Indeed, numerical evaluations of the 
function b(z) at the end of J9j and in lfl4l suggests to us that a moderate anisotropy would 
only slightly modify an optimal geometry found in the isotropic limit. We leave, therefore, the 
full investigation of anisotropic effects for future research. Finally, the quantum efficiency of 
exciton generation is taken to be unity in both donor and acceptor regions: <p, ■ = 1 . 

Given geometrical parameters and refraction indexes at A ( , the function g(z) can easily be 
computed with standard transfer matrix methods ll35ll40l . Assuming a unit-amplitude wave 
exp (— in ext koz) coming from z = °°, one finds the wave a ; exp(— in^koz) + £>,exp (m^fcoz) in 
region ; and 

Re(n (> l ■ ■ ■ . 2 



(23) 
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The numerical evaluation of integrals of the type which appears in (fT~5b and (l20l may present 
some difficulties when the denominators of the integrand approach zero. To avoid this problem, 
it was recommended to deform the integration contour away from these poles l26l . This is 
especially interesting when a large number of evaluations are required and when the poles are 
numerous. In our case, we checked that it was unnecessary, as the number of poles was small 
and that they were in practice sufficiently far from the real axis. More troublesome, numerically, 
can be the square root singularities at k e ,\ = and k a -\ = (both expressions being identical if 
e x = e z .) These are simply removed by changing the integration variable from k^ to either k e ^\ 
or k ,\ . Furthermore, it is useful to adopt the wavenumber unit such that Icq = 1 . Regarding the 
integration of (O, we discretized each photoactive region with a Chebychev grid BP . In this 
way EQE SC was accurately computed with significantly fewer discretization points, and hence 
fewer evaluations of b(z), than with standard finite differences. 

3.1. Parallel excitons 

Based on the geometry considered in the introduction, we first look for a structure that is fa- 
vorable to parallel excitons. It was noted that a low refractive index next to silver can bring 
b\\ (z) to a minimum near the back electrode, but that b\\ (z) rapidly increases further away from 
it. It is therefore desirable to extend the range of low b\\ values away from the back electrode 
so as to combine long lifetime with good sunlight harvesting. This is achieved with the con- 
figurations described in Table Q] One configuration uses ITO as the transparent electrode, the 
other one uses Ag. The parameters were determined to maximize EQE SC for q = 1 and a bulk 
diffusion length of 5 nm. Figs. [3^ and b show the corresponding functions b\\(z), b±(z), and 
g(z). One may see that, indeed, low b\\ values are maintained over a longer distance from the 
back electrode than in Drexhage's configuration. Moreover, and contrary to intuition, Fig. [3^ 
gives an example where both b« and b± are small at the same place. Solving (|5]l and (|6]l with 
these parameters and various bulk diffusion lengths, L, we plot the EQE SC as a function of L 
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Fig. 3. Decay rates (q = 1), gain functions, and external quantum efficiencies for the param- 
eters given in Table [T] Top: ITO transparent electrode. Bottom: Ag transparent electrode, 
iso: isotropic case where b = + %b±. 



Table 1. Two configurations optimized for parallel excitons. A sunlight absorption 
length of 70 nm is assumed in both photoactive materials. ;iito = 1.76 + 0.08/, 
«A g (900nm) = 0.04 + 6.37/, K Ag (750nm) = 0.03 + 5.19/. 
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in Figs. |3j; and d. This allows one to compare the merit of non-fluorescent {q = 0) and purely 
fluorescent (q = 1) materials with parallel, perpendicular or randomly oriented excitons. The 
theoretical advantage of q = 1 appears clearly, even in the isotropic case. This is particularly 
true when the diffusion length is small compared to the thickness of the active region. Most 
of the q — 1 curves have almost reached their asymptotic value for L = 5 nm while twice that 
diffusion length is necessary for q = 0. Hence, the effective diffusion length is approximately 
doubled in these examples thanks to fluorescence. 

All the curves EQE SC (L) displayed are very well approximated by ([TO with suitable values of 
b estimated from Figs[3^,b. The ultimate EQE SC , as L becomes large, is given by the absorption 
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Fig. 4. Absorption efficiency T}a computed for the devices described in Table Q]as a func- 
tion of sun wavelength. A fixed ITO refractive index nrro = 1 -76 + 0.08/ and an absorption 
length of 70 nm in the active region are assumed over the whole spectral range. Ag refrac- 
tive index is taken from |34|. 

efficiency, v\a- At the wavelength considered, the use of an Ag transparent electrode appears 
advantageous compared to ITO, as it allows one to obtain large values of v\a- If one varies 
Aj over some tens of nm, however, one witnesses a sharper drop of 77^ with Ag than with ITO. 
Hence, a silver transparent electrode makes light injection more resonant than an ITO electrode. 
This is illustrated in Fig. [4] 

3.2. Perpendicular excitons 



Table 2. Optimized configuration for perpendicular excitons. Absorption length: 70 
nm is assumed in both photoactive materials. rarro = 1-76+0.08/, nAg(900nm) =0.04 + 
6.37/, n Ag (750nm) = 0.03 + 5. 19/. 
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We now turn to the transport of perpendicular excitons. Considering the formula for b±(z) 
in (T5i , we note that large absolute values of Rq and R'{ are required in order to significantly alter 
the decay rate. A way to favor such a situation is to raise the refraction index in the photoactive 
region. On the other hand, a large index contrast may induce a strong reflection of the sunlight 
at the D/EBL interface, which is detrimental to the gain function g(z). A compromise between 
the two effects is given in Table|2] As with parallel excitons, we find that a metallic transparent 
electrode can give rise to a higher short-circuit current than with ITO and we thus focus on 
a silver electrode. The gain and loss profiles corresponding to Table [2] are shown in Fig. 
Remarkably small decay rates are obtained in the photoactive region, resulting in a dramatic 
improvement of EQE SC when passing from q = to q = 1 . 

The very small decay rates shown in Fig.|5]and in Fig. [3^ motivate us to try and find a simple 
estimate of the minimum achievable b± in a solar cell. To this end, let us assume a high-index 
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Fig. 5. Gain and decay profiles computed with Tableland q = 1. HBL: hole-blocking 
layer, A: acceptor, D: donor, EBL: electron-blocking layer. 



photoactive material (ri\) between low-index blocking layers («2). Restricting our attention to 



isotropic media, we have £,• x = £,• - = rq and k ez j = k„ 



k z j. In particular, 



k 7 . 



K :.2 



(24) 



On account of the thinness of the photoactive region, we neglect the phase factors in Rq and R P y 
and approximate them by R% and R^ , respectively. Furthermore, we assume that the electrodes 
are sufficiently remote for their influence to be negligible. This point is partly justified for values 
of fen in (TT3T > such that k Z: 2 is imaginary but is not otherwise rigorous. Thus, we write 



and the expression for b± becomes 
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Using d24l and noting that jc 3 / (a — x 2 
can be evaluated analytically, yielding 
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(25) 



(26) 



(27) 
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the above integral 



n 2 

q\ — 



(29) 



In particular, for unit fluorescence quantum yield, we have the approximate formula b± w 
(i2/«i) 5 - As shown in Fig. [6] this simple expression provides a fairly good approximation for 
the whole range of index contrast allowed by 1 .3 < «2 < "l < 2.8. With n\ = 2.8 and 112 = 1 .3, 
the minimal achievable decay rate is estimated as £>x, m in ~ 0.022, to be compared with the 
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Fig. 6. b^(z) in a structure with refractive indexes ;2A g /«2/ ra l / n 2/ n Ag/ n 3 an d thick- 
nesses (nm): 140/29/295/6. n 2 = 1.3, n } = 2.8. Cases a to e: n x = 2.8,2.2, 1.7, 1.5, 1.4. 
Wavelength: 900 nm. Dashed line: numerical. Full line: («2/"i) 5 - 



actual value, 0.027, in Fig. [6] Moreover, ((29} remains a good estimate as long as the distance to 
the nearest silver electrode exceeds 100 nm. Below that distance, energy transfer to the silver 
electrodes, notably through surface plasmon polaritons H42I . becomes important. A remarkable 
feature of the approximation ((29) is that it is independent of the exciton wavelength and of 
geometrical parameters. 

4. Conclusion and perspectives 

Next to the usual material characteristics such as absorption spectrum, LUMO and HOMO, ex- 
citon binding energy, lifetimes and mobilities, this study suggests to also consider fluorescence 
quantum yield, q, in the quest for good photovoltaic molecules. Molecules with a high quan- 
tum yield are amenable to fluorescence management, by which the exciton diffusion length can 
be considerably increased. The approximate formula (TPH indicates that, for thin photoactive 
regions, 

EQE sc = 77 a (A. v )tj d (A, ? ), (30) 

which brings out the necessity to jointly consider the light harvesting capability (rj^) and the 
transport property of a solar cell (?]£>). Both depend on the geometry, but fluorescence affects 
tj d only. 

Our study contradicts several expectations derived from Drexhage's configuration. To start 
with, the strongest fluorescence inhibition was found for perpendicular excitons, rather than 
for parallel ones. An approximate rule is that perpendicular exciton lifetime can be increased 
up to a factor (ni/«2) 5 , where n\ is the refraction index in the photoactive material and ri2 is 
the index in the neighboring layers. This amounts to increasing the diffusion length by a factor 
[nxjn^y'l 1 with respect to the bulk value. This formula holds well if the HBL and EBL are 
sufficient thick for dissipative energy transfer to be negligible and if the active region is thin 
compared to the fluorescence wavelength. Secondly, by an appropriate choice of the optical 
and geometrical parameters, one can have reduced fluorescent decay rate for perpendicular and 
parallel excitons at the same places. Hence, improvement of device performance is theoretically 
possible for q — 1 also with random exciton orientation. 

The ultimate EQE SC is 77,4. In this respect, a metallic transparent electrode such as Ag was 



found to yield larger values than a dielectric one (at optical frequencies), such as ITO. Moreover, 
in the example given in Table [T] an Ag thickness of 13 nm was found to be preferable than the 
more usual 10 nm. Indeed, the resulting fluorescence inhibition and improved absorption in the 
photoactive region at X s = 750 nm more than compensate absorption losses in the transparent 
electrode incurred by passing from 10 to 13 nm thickness. At the same time T}a(X s ) shows a 
sharper spectral resonance with Ag than with ITO, see Fig. |4] For simplicity, Fig. |4] was drawn 
by assuming a fixed absorption length over the whole frequency range considered. In practice, 
the most judicious choice between ITO and a metallic transparent electrode should involve 
matching the microcavity response with the molecular absorption spectrum. 

The optimal geometries presented in this paper are slightly unusual, particularly regarding the 
large thicknesses of the HBL and EBL. Such thicknesses are not uncommon in OLED design. 
Presently, they are necessary to make fluorescent excitons immune to dissipative energy transfer 
in the electrodes. 

Our work also shows that, for a badly designed geometry, a purely fluorescent material is 
less efficient than a non-fluorescent material with the same bulk lifetime. This can easily be 
seen by considering Fig. [3] If the gain material was located between 150 and 200 nm in that 
example, then both b± and bn would be larger than one. In that case, obviously, rjo(X,q^ 1) < 
r\o(X : q = 0). This may be surprising, as Shockley and Queisser showed that a large quantum 
yield was necessary to maximize the solar cell efficiency ||43ll44ll . However, the Shockley- 
Queisser limit was derived for the much thicker semiconductor devices, where the microcavity 
effects discussed here are negligible and where exciton transport is not an issue. 

In this work, we have only made a partial exploration of the parameter space. There is there- 
fore still room for improvement or adaptation to technical realities. In the frame of the present 
theory, the number of configurations that can be studied exceeds by far the examples that we 
have discussed. We have limited our parameter investigation to isotropic materials, to donor 
and acceptor molecules having unit quantum yield, identical diffusion lengths, and identical 
orientation. These constraints can obviously be lifted with the mathematical model presented 
here. Besides, other sunlight injection strategies could possibly be devised that combine large 
value of 7]a with large 7]d- 

Finally, although we focused on large-g molecules, the present theory should also be useful 
to model materials with moderate values of q. Indeed, the functions b± and bn can still vary 
significantly from unity in that case. Models like (fJJ then appear to be more sensible than ((TJ, 
while still simple to implement. In particular, important experimental estimations of diffusion 
lengths based on the assumption of spatially uniform exciton lifetimes may have to be revised 
for some molecules. 
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A. Appendix: Derivation of (Il5l and d20t 

In this appendix, we derive formula (fTBl and (l20l for a uniaxial multilayer. As a preample, let 
us note that if the purely fluorescent part of the exciton decay rate can be written as 



where F, is the bulk value and y(z) is the correction due to the boundaries, then the total decay 
rate is 



r r (i + y(z)), 



(31) 




(32) 



Therefore, 

b(z) = l+qy(z). (33) 

Considering time-harmonic electromagnetic fields, Ampere's and Faraday's law can respec- 
tively be written as 

fcgeE = -i£»/iJ + io)VxB, (34) 
z<»B = V x E, (35) 

where Icq = ^/EqI-ICO and £ is the permittivity tensor, of the form given in (fTZt . Note that ( l33T > 
automatically makes B divergence-free and that V • D = p follows from d34i l and the continuity 
equation V • J icop = 0. We find it useful, as in [45 ,46 1, to use the z components of E and B as 
electromagnetic potentials. We shall decompose any vector Z into its z component and its (x,y) 
projection as 

Z = Z||+zZ-. (36) 
After some manipulations, we find that the plane components of ( f34b and d35l ) can be written as 

d 2 ,\ „ dE, 



+ sM En = — i'couj + Vii it - icoz x V«B 7 , (37) 
az J oz 

( d 2 ,\ , . d . s dB, 7 

I -j^+BxK j (i«)B||) = -icofj.-^ (zx J||) +ifi)V||-^-*oe r zx V||£,, (38) 



while 



£,-^2 + fcVfl + £ T £-^J £, = -i(a\ie x z ■ \i + ^ VV ■ J ) , 



(^2 + V f + e -^o) = -Mz • (V X J) 



(40) 



Equations (T39b and (l40b make it apparent that E z is the electromagnetic potential for extraordi- 
nary waves and that B z is the potential for ordinary waves. Note that the left-hand-side of (T39b 
can be transformed into the Helmholtz equation by rescaling the coordinates. Therefore, the 
Green's functions for Eqs. ((39) and (l40l in free space are, respectively [46], 



go = — -====-. (41) 
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These functions have the following representation, due to Sommerfeld 



- f^M^Wibdih,, g0 = ^ r J 4^e*°M hdkh (42) 
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where p = yx 2 + y 2 and k e . z , k„. z are given by the dispersion relations (Q~3) and (Il4l in the 
medium being considered. 

With J = jo<S(x), we may now write the solution of (|39l and d40l > as 

E z = -i(o^e x (^z-jo + ^y z io-^ge+E', B z = iixz-{j xV)g o +B', (43) 

where E' and B' are non singular solutions of the homogenous differential problem that ensure 
that boundary conditions at finite distance are satisfied. 



A. 1. Perpendicular dipole 

Let us assume first that jo = joz. We may then take B z = and, combining (f42t and d43l l. write 

E- as 

E z = -P^g f JL/ (jfc, |P ) f^W +C e^+De-^A dk h (44) 



where C and D are arbitrary functions of h\ . Let the nearest interfaces be at a distance l\ above 
and a distance £o below the dipole; these are z\— z and z — zo in ( Tl3T > and ( l20l i. The boundary 
conditions at these interfaces translate into reflection coefficients for the waves that compose 
the integral in (l44l . For a given k\\, a wave with amplitude (1 +C) travels upwards from the 
origin. It is reflected with a coefficient R p at z = l\, to be computed according to (fT8T l. Hence, 
D should be 

£> = R\e m '^ (1 + C) = i?f (1 + C) . (45) 

By the same token, considering waves traveling downwards from the origin to the interface at 
z = —£q, the amplitude C is given by 

C = R p e 2iKJo (l + D) = R p {l+D). (46) 

Solving d43T l and d46l i. we thus find that 

*SM+f& R P R P +R P 
r = 010 D = 1 — l - (47) 

l-R p R p ' l-R p R p ' K } 



and, hence, that 



With this solution, Eq. ( f37l > for E|| becomes 
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The solution is, simply, 
E =_^Lv d r k \\Mk\\P) (j ke . r \z\ , ^i+*o.fe , W+^-a^ * 
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where the last term in brackets ensures the continuity of E|| at z — 0. The purely fluorescent 
decay rate is given by the rate of work done by the dipole divided by the photon energy: 

-2Re{E-jg} = V£x\jo\ 2 Rc \ r k \ ( 1 Iftm+Zj + K 
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The bulk value is 



Tr = E^[ £l ^± Ah = ^^ell\ (52) 
2%%em Jo ke iZ " 3>nh 

Using this expression, we may rewrite the right hand side of ([5Tb as 
/ 3e! /2 f r k l 

from which (JT3J follows. 




A.2. Parallel dipole 

Let now jo = j'ox. This time, both ordinary and extraordinary waves are generated. Following 
the same reasoning as for the perpendicular dipole, one obtains 

-conk ,? 2 rMk\\p) ( ikeM , K&l-K , KK-K 

Z I /til Lift 

(54) 

B = u i ^u^^ , gikoM »o»i -r»o e ,i ,c + fVM^l e -iko# ) k dk . (55) 
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With these solutions, we may write the equation for E|| . Noting that 

5( X ) = ^M/' J (k«p)h\dk h (56) 
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we have 



' e ^o )E|| =^^x5(z) / y (*||P)*||d*|| 
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^MzxV,^ r*M f^W + S^^ + S^^ih^,. (57) 
4tt 11 dyJo k . \ 1-RqR\ l--RS R i / 

Given that (2/A:) -1 exp(//:|z|) is the solution of +k 2 f = <5(z), the solution of (|57] > is 
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which simplifies into 
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We may now compute the fluorescent decay rate (at z = 0,p = 0): 
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The bulk value, T, is obtained when all reflection coefficients vanish: 
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Factorizing this expression out of (l60t . we obtain 
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and < f20b follows. 



